Symmetry energy of deformed neutron-rich nuclei 
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The symmetry energy, the neutron pressure and the asymmetric compressibility of deformed 
neutron-rich even-even nuclei are calculated on the examples of Kr and Sm isotopes within the 
coherent density fluctuation model using the symmetry energy as a function of density within the 
Brueckner energy-density functional. The correlation between the thickness of the neutron skin and 
the characteristics related with the density dependence of the nuclear symmetry energy is investi- 
gated for isotopic chains of these nuclei in the framework of the self-consistent Skyrme-Hartree-Fock 
plus BCS method. Results for an extended chain of Pb isotopes are also presented. A remarkable 
difference is found in the trend followed by the different isotopic chains: the studied correlations 
reveal a smoother behavior in the Pb case than in the other cases. We also notice that the neutron 
skin thickness obtained for 
data. 



Pb with SLy4 force is found to be in a good agreement with recent 
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I. INTRODUCTION 

The study of the nuclear matter symmetry energy that 
essentially characterizes the isospin-dependent part of the 
equation of state (EOS) of asymmetric nuclear matter 
(ANM) is currently an exciting topic of research in nu- 
clear physics [U-iI|. In fact, applications of ANM are 
broad, ranging from the structure of rare isotopes [H [f| 
to the properties of neutron stars @, H| and the dynami- 
cal process of nuclear reactions @. The transition from 
ANM to finite nuclei is a natural and important way to 
learn more about the nuclear symmetry energy which is 
poorly constrained by experimental data on ground-state 
nuclear properties. 

The ground states of atomic nuclei are characterized 
by different equilibrium configurations related to corre- 
sponding geometrical shapes. The study of the latter, as 
well as the transition regions between them, has been a 
subject of a large number of theoretical and experimen- 
tal studies (for a review, see, for example, Ref. [l(J and 
references therein) . The position of the neutron drip line 
is closely related to the neutron excess and the deforma- 
tion in nuclei. Deformed nuclei are expected in several 
regions near the neutron drip line [TTI.[r3|. In some cases, 
the deformation energy can impact their existence. For 
instance, it has been predicted that there exist particle- 
bound even-even nuclei that have, at the same time, neg- 
ative two-neutron separation energies caused by shape 



coexistence effects In fact, the nuclear deformation 
increases the surface area, thus leading to a larger sur- 
face symmetry energy in a neutron-rich nucleus with a 
deformed shape. Conversely, the precise determination 
of the surface symmetry energy is important to describe 
the deformability of neutron-rich systems and also to val- 
idate theoretical extrapolations. 

The theoretical treatment of the ground-state prop- 



erties of deformed nuclei is usually made in the frame- 
work of mean-field approaches, for instance, the non- 
relativistic self-consistent Hartree-Fock+BC S appro ach 
with density-dependent Skyrme interaction [13l4l7| . or 
the Hartree-Fock-Bogoliubov (HFB) approximation that 
uses, in particular, the Gogny force containing a finite- 
range interaction [l^, [l9| . Also the relativistic mean- field 
(RMF) models with different types of parameter sets 
have been widely used. For instance, the RMF theory 
with parameter set FSUGold has been recently employed 
to study the binding energies, quadrupole deformations, 
charge radii and neutron skins of rare-earth even-even 
nuclei ranged from Z = 58 to Z = 70 and some de- 
formed nuclei (Nd, Sm, Gd, Dy) in other regions [20| . 
This new parameter set that includes the nonlinear cou- 
pling between the isoscalar and isovector mesons is shown 
to reproduce successfully the ground-state properties of 
deformed nuclei. Also, the same parameter set FSUG- 
old can successfully reflect the shell effect of the neutron 
magic number N = 82 thus leading to a good descrip- 
tion of t he g round-state properties of Sn, Te, Xe, and Ba 
isotopes |21j |. 

Nowadays, the experimental information about the 
symmetry energy is fairly limited. The need to have in- 
formation for this quantity in finite nuclei, even theoret- 
ically obtained, is a major issue because it allows one to 
constrain the bulk and surface properties of the nuclear 
energy-density functionals (EDFs) quite effectively. For 
example, the traditional Skyrme EDF based on the lep- 
todermous expansion of the smooth nuclear energy was 
used in Ref. [22J, where the strongly correlated symme- 
try and surface symmetry terms in this expansion are 
resolved by considering data on deformed neutron-rich 
nuclei in which the surface symmetry term is amplified. 
By starting from microscopic nucleon self-energies in nu- 
clear matter, the authors of Ref. [23l have made a com- 
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prehensive study of deformation properties of relativistic 
nuclear EDFs for a set of 64 axially deformed nuclei in 
the mass regions A «150-180 and A «230-250. 

In our recent work [24| the Brueckner EDF for infinite 
nuclear matter was applied to calculate nuclear quanti- 
ties of medium-heavy and heavy Ni, Sn, and Pb nuclei 
that include surface effects, namely the nuclear symme- 
try energy s, the neutron pressure po, and the asym- 
metric compressibility AK. For this purpose, a theo- 
retical approach that combines the deformed HF+BCS 
method with Skyrme-type density-dependent effective in- 
teractions 1 13] and the coherent density fluctuation model 
(CDFM) [Mill was used. We would like to note the 
capability of the CDFM to be applied as an alternative 
way to make a transition from the properties of nuclear 
matter to the properties of finite nuclei. We have found 
that there exists an approximate linear correlation be- 
tween the neutron skin thickness AR of even-even nuclei 
from the Ni (A = 74 - 84), Sn (A = 124 - 152), and 
Pb (A = 206 — 214) isotopic chains and their nuclear 
symmetry energies. A similar linear correlation between 
AR and po was also found to exist, while the relation 
between AR and AK turned out to be less pronounced. 
The kinks displayed by Ni and Sn isotopes and the lack 
of such kink in the Pb chain considered [24{ were shown 
to be mainly due to the shell structure of these exotic 
nuclei but they deserve further analysis within the used 
theoretical approach. 

Another interesting question is to explore how the nu- 
clear symmetry energy changes in the presence of defor- 
mation and correlates with the neutron skin thickness 
within a given isotopic chain. In Ref. [llj], the effects of 
deformation on the skin formation were studied in Kr iso- 
topes that are well-deformed nuclei. It has been shown 
from the analysis on 98 > 100 Kr nuclei that although the 
profiles of the proton and neutron densities, as well as 
the spatial extensions change with the direction in both 
oblate and prolate shapes, the neutron skin thickness re- 
mains almost equal along the different directions perpen- 
dicular to the surface. Thus, a very weak dependence of 
the neutron skin formation on the character of the defor- 
mation was found [l6j . 

In the present work an investigation of possible re- 
lation between the neutron skin thickness and the ba- 
sic nuclear matter properties in deformed finite nuclei, 
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such as the symmetry energy at the saturation point, 
symmetry pressure, and asymmetric compressibility, is 
carried out for chains of deformed neutron-rich even- 
even Kr (A = 82 — 96) (including, as well, the case 
of some extreme neutron-rich nuclei up to 120 Kr) and 
Sm (A = 140 — 156) isotopes, following the theoretical 
method of Ref. [24] . We also present for comparison re- 
sults for an extended chain of Pb (A = 202 — 214) iso- 
topes. This is motivated by the si gnifi cant interest (in 
both experiment [2714291 ] and theory j30l - [33l |) to study the 
neutron distribution and rms radius in 208 Pb, aiming at 
precise determinations of the neutron skin in this nucleus. 
In addition to the interest that this study may have by 
itself as well as in combination with the previous calcu- 
lations of Ref. 24j, we give some numerical arguments 
in proof of the existence of kinks in Ni and Sn isotopic 
chains that are not present in the Pb chain. The kinks 
are produced because of the sensitivity of the symmetry 
energy and neutron pressure to the shell structure (see, 
for instance, the discussion in Refs. 0, HU, HH|). 

The paper is organized as follows. In Sec. II we 
present a brief description of the theoretical formalism 
(definitions of ANM properties, CDFM basic expressions, 
Brueckner energy-density functional, Hartree-Fock+BCS 
densities) used to unveil a possible correlation between 
the neutron skin thickness and the nuclear matter char- 
acteristics of the considered isotopic chains. Section III 
contains our results with a discussion on the obtained re- 
lationships and on the presence of kinks. The concluding 
remarks are drawn in Sec. IV. 



II. THEORETICAL FRAMEWORK 

We study in the present work the symmetry energy 
s(p) and related quantities of finite deformed nuclei on 
the basis of the corresponding definitions for ANM. The 
quantity S ANM (p) ( which refers to the infinite system and 
therefore neglects surface effects, is related to the second 
derivative of the energy per particle E(p, 5) using its Tay- 
lor series expansion in terms of the isospin asymmetry 
S = (p n — Pp)/p {p, Pn and p p being the baryon, neutron 
and proton densities, respectively) (see, e.g., [24|, HH, H3 ) : 



s ANM (p) = - d2E ( p > 6 "> 
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p ANM , . AK ANM , 
= a 4 + — [p - pa) + (p - po) + 

s=o Po i! Vo 



(1) 



In Eq. (JlJ the parameter 0,4 is the symmetry energy at 
equilibrium (p = po). In ANM the pressure p ANM and 
the curvature AK 



ANM 



ANM _ 2 

Po — Po 



are: 

ds ANM 



(P) 



dp 



(2) 



P=Pa 



AK ANM = Qf) 2 



d 2 S ANM (p) 



dp 2 



(3) 
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The "slope" parameter L 



ANM 



is defined as 



- ANM 



ANM 



Po 



(4) 



In general, the predictions for the symmetry energy vary 
quite substantially: e.g., 04 = s(p ) = 28 — 38 MeV while 
an empirical value of 04 « 29 MeV has been extracted 
from finite nuclei by fitting the ground-state energies us- 
ing the generalized Weizsacker mass formula (see, e.g., 
Ref. [HI). By using the experimental pygmy strength, an 
average value of 0,4 = 32.0 ± 1.8 MeV was obtained from 
the 130 ' 132 Sn analysis [3(|, which is within the acceptable 
range of values of 04 to be around 32.5 MeV coming from 
various experiments using different experimental probes 
(for a recent status, see, for example, Ref. [5(| and refer- 
ences therein). 

In Ref. [24j | we calculated the symmetry energy, the 
pressure and slope, as well as the curvature for finite nu- 
clei applying the coherent density fluctuation model (sug- 
gested and developed in Refs. [la El]). In the CDFM the 
one-body density matrix p(r, r') of the nucleus is written 
as a coherent superposition of the one-body density ma- 
trices p x (r,r') for spherical "pieces" of nuclear matter 
called "fluctons" with densities p x (r) = Po(x)Q(x — |r|), 
p (x) = 3A/4irx 3 : 



p(r,r')= / dx\f(x)\ 2 p x (r,r') 
Jo 



(5) 



with 



p x {r,r') = 3p (x) 
x 6 I x- 



^(M^lr-r'l) 
(k F {x)\v-v'\) 
|r + r' 



where j\ is the first-order spherical Bessel function, 
'3tt 2 ^ 1/3 



k F{x) = ( — p (x) 



a 
x 



with 



1/3 



1.52^4 1 / 3 



(6) 



(7) 



(8) 



is the Fermi momentum of the nucleons in the " fiucton" 
with a radius x. In Eq. §5§ \f(x)\ 2 is the weight func- 
tion that in the case of monotonically decreasing local 
densities (dp{r)/dr < 0) can be obtained using a known 
density distribution for a given nucleus: 



1 



dp(r) 



Po(x) dr 



(9) 



with the normalization J °° dx\f(x)\ 2 = 1. 

The main assumption of the CDFM is that properties 
of finite nuclei can be calculated using the correspond- 
ing ones for nuclear matter, folding them with the weight 



function |/(a;)| 2 . Along this line, in the CDFM the sym- 
metry energy for finite nuclei and related quantities are 
assumed to be infinite superpositions of the correspond- 
ing ANM quantities weighted by |/(cc)| 2 : 



dx\f(x)\ 2 s ANM {x) 



Po = I dx\f(x)\ 2 p ANM (x), 
10 



AK = 



dx\f(x)\ 2 AK ANM (x). 



ANM 



(10) 



(11) 



(12) 



The explicit forms of the ANM quantities s 
p ANM (x), and AK ANM (x) in Eqs. (HUJ), @D, and U4> 
are defined below. They have to be determined within a 
chosen method for description of the ANM characteris- 
tics. In the present work, as well as in Ref. [24j |. consid- 
ering the pieces of nuclear matter with density po(x), we 
use for the matrix element V(x) of the nuclear Hamilto- 
nian the corresponding ANM energy from the method of 
Brueckner et al. [HI H2] : 



V(x) =AV {x)+Vc-V co , 



(13) 



where 



V (x) 



37.53[(1 + S) 5/3 + (1 - 5) 5/3 }p 2 /3 (x) 
hpo(x) + b 2 Po /3 (x) + b 3 pl /3 (x) 
S 2 [b iPa (x) + b 5 p 4 /3 (x) + b 6P 5 /3 (x)] (14) 



with 



h = -741.28, b 2 = 1179.89, 63 = -467.54, 

bi = 148.26, 65 = 372.84, b 6 = -769.57. (15) 

In Eq. (1131) Vo(x) is the energy per particle in nu- 
clear matter (in MeV) accounting for the neutron-proton 
asymmetry, Vc is the Coulomb energy of protons in a 
fiucton, 



V c 



3 Z 2 e 2 
5 x 



and Vco is t ne Coulomb exchange energy: 
V CO = 0.7386Ze 2 (3Z/47ra 3 ) 1/3 . 



(16) 



(17) 



Thus, using the Brueckner theory, the symmetry energy 
x) and the related quantities for ANM with den- 



ANM, 



sity po(x) (the coefficient 04 in Eq. (fT])) have the forms: 

s ANM (x) = A1.7p 2 /3 (x) + b 4 p (x) 

+ b 5 p 4 /3 (x) + b 6 pl /3 (x), (18) 



p ANM (x) = 27.8p 5 /3 (x)+b iP 2 (x) 

1 7 /3( \ , 5 l 8/3/ ^ 



(19) 
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and 



AK ANM (x) = -83 Ap 2 J\x) + 4b 5P 4 /3 (x) 
+ I0b 6 p 5 /3 (x). 



(20) 



In our method (see also Eqs. (H]), (H]), and (J2UJ) arc 
used to calculate the corresponding quantities in finite 
nuclei s, p , and AK from Eqs. (10]), (11]), and (1T21) . 
respectively. We note that in the limit case when p(r) = 
poQ(R-r) and |/(a;)| 2 becomes a S function [see Eq. ©], 
Eq. (TO]) reduces to s ANM {p ) = a 4 . 

In our work we use the proton and neutron densi- 
ties obtained from self-consistent deformed Hartree-Fock 
calculations with density-dependent Skyrme interactions 
[l3j and pairing correlations. Pairing between like nu- 
cleons is included by solving the BCS equations at each 
iteration with a fixed pairing strength that reproduces 
the odd-even experimental mass differences [43j . 

The spin-independent proton and neutron densities are 
given by [H 



p(R) = p(r,z) 



(21) 



where r and z are the cylindrical coordinates of R, vf 
are the occupation probabilities resulting from the BCS 
equations and pi are the single-particle densities 



Pi (R) = Pi (r, z) = |$+(r, z)\ 2 + |$r( r , z)\ 2 (22) 



with 



'2tt 



x Y, fe,±i/2VQ<M^W (23) 

OL 

and a = {n r , n z , A, £}. In (23]) the functions if>^ (r) and 
■0„ z (z) are expressed by Laguerre and Hermite polyno- 
mials: 



0) = 



(n r + A)! 



P ± V2r, A / 2 e-^ 2 Lt(v), (24) 



^«.(*) = J r L ^ l' 2 e- f/2 H n M (25) 
y V7r2"-n 2 ! 

with 

(3 Z = (mujz/h) 1 / 2 , ± = (mu J± /h) 1 / 2 , 

C = z(3 z , r/ = r 2 /3 2 . (26) 

The normalization of the densities is: 



p(R)dR = X (27) 
with X = Z, N for protons and neutrons, respectively. 



The multipole decomposition of the density can be 
written in terms of even A multipole components as 

mm 



p(r,z) =Y P\ {R)P\ (cos 9). 



(28) 



In the calculations, for the density distribution p(r) 
needed to obtain the weight function |/(ir)| 2 [Eq. (9|] 
we use the monopole term po(R) in the expansion (28]). 

The neutron skin thickness is usually estimated as the 
difference of the rms radii of neutrons and protons: 



AR 



--< rl > 1 / 2 



< r 2 p > 1 / 2 



(29) 



In our calculations the following Skyrme force 
parametrizations are used: SLy4 [Hj], Sk3 SGII (47j . 
and LNS [48|. These are among the most extensively 
used Skyrme forces that work successfully for describing 
finite nuclei properties. Although it is well known that 
SGII and LNS interactions do not predict accurate 
binding energies in finite nuclei, we have included them 
in our work because they are representative examples of 
Skyrme interactions and because we are not concerned 
here with the absolute values of binding energies, but 
rather with the isotopic evolution of relative differences 
of various magnitudes and energy derivatives. 



III. RESULTS AND DISCUSSION 

We start our analysis by searching for the role of de- 
formation on the rms radii. It is illustrated by showing 
the differences between the proton (neutron) radii and 
the corresponding proton (neutron) radius of semi-magic 
86 Kr (N = 50) which is taken as a reference nucleus. The 
results are given in Fig. Q] [(a) and (b), respectively]. In 
addition, in Fig. QJc) our results with SLy4 force for the 
squared charge radii differences in Kr isotopes are shown 
and they are compared with the experimental data from 
Ref. Es]. Results of this type have been presented in 
Ref. [161 ]. where the charge radii differences in Sn iso- 
topes obtained from SLy4, SGII, and Sk3 Skyrme forces 
are compared with the experimental data taking the ra- 
dius of 120 Sn as the reference. In principle, these dif- 
ferences have been found to be very sensitive probes of 
nuclear shape transitions and it is worth studying glob- 
ally a wide region of neutron-rich exotic nuclei and dis- 
cussing the similarities and differences among the various 
isotopic chains [50j . 

In this analysis the even-even 86 Kr isotope turns out to 
be spherical, while for Kr isotopes beyond the semi-magic 
86 Kr (N = 50) nucleus a deformation (oblate or prolate) 
takes place. We note that for the considered lightest iso- 
topes with A = 82, 84 both shapes produce results that 
are indistinguishable and, therefore, only the prolate so- 
lutions are presented. A smooth increase of radii differ- 
ences relative to 86 Kr with increasing neutron number 
can be observed in all three panels of Fig. [T] The charge 
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and proton radii grow following the growth of the neutron 
radius with increasing neutron number. These increases 
are only slightly dependent on whether an oblate or a 
prolate shape is considered. A similar behavior has been 
obtained from Gogny-DlS-HFB calculations performed 
in Ref. (sjj. A satisfactory agreement with the experi- 
mental isotope shifts is observed in Fig. HJc) that pro- 
vides a good starting point to study further quantities 
such as the symmetry energy and related characteristics 
of deformed nuclei within our theoretical method. 

Next, an illustration of a possible correlation of the 
neutron-skin thickness Ai? with the s and po parameters 
extracted from the density dependence of the symmetry 
energy around the saturation density for the Kr isotopic 
chain is given in Fig. [5J The symmetry energy and the 
pressure are calculated within the CDFM according to 
Eqs. (fT0|) and (fTTj) by using the weight functions (|9|) cal- 
culated from the self-consistent densities in Eq. (|2"Tj) . The 
differences between the neutron and proton rms radii of 
these isotopes [Eq. (|29|) ] are obtained from HF+BCS cal- 
culations using four different Skyrme forces, SLy4, SGII, 
Sk3, and LNS. It can be seen from Fig. [5] that there ex- 
ists an approximate linear correlation between Ai? and 
s for the even-even Kr isotopes with A = 82 — 96. Sim- 
ilarly to the behavior of Ai? vs s dependence for the 
cases of Ni and Sn isotopes [2~jj . we observe a smooth 
growth of the symmetry energy up to the semi-magic nu- 
cleus 86 Kr (N = 50) and then a linear decrease of s while 
the neutron-skin thickness of the isotopes increases. This 
linear tendency expressed for Kr isotopes with A > 86 is 
similar for the cases of both oblate and prolate deformed 
shapes. We note that all Skyrme parametrizations used 
in the calculations reveal similar behavior; in particular, 
the average slope of Ai? for various forces is almost the 
same. 

In addition, one can see from Fig. [2] a stronger devia- 
tion between the results for oblate and prolate shape of 
Kr isotopes in the case of SGII parametrization when dis- 
playing the correlation between Ai? and s. This is valid 
also for the correlation between Ai? and po , where more 
distinguishable results for both types of deformation are 
present. The neutron skin thickness Ai? for Kr isotopes 
correlates with po almost linearly, as in the symmetry- 
energy case, with an inflexion point transition at the 
semi-magic 86 Kr nucleus. In addition, one can see also 
from Fig. [5] that the calculated values for pq are smaller 
in the case of LNS and SLy4 forces than for the other 
two Skyrme parameter sets. In general, we would like to 
note that the behavior of deformed Kr isotopes shown in 
Fig. [5] is comparable with the one found for the spher- 
ical Ni and Sn isotopes having a magic proton number 
that we discussed in Ref. [24[ . The small differences just 
indicate that stability patterns are less regular within 
isotopic chains with a non-magic proton number. 

For more complete study, we also consider in our work 
the extremely neutron-rich Kr isotopes (A = 96 — 120). 
The results for the symmetry energy s as a function 
of the mass number A for the whole Kr isotopic chain 



(A = 82 — 120) are presented in Fig. [3] We observe peaks 
of the symmetry energy at specific Kr isotopes, namely 
at semi-magic 86 Kr (N = 50) and 118 Kr (N = 82) nuclei. 
In addition, a flat area is found surrounded by transi- 
tional regions A = 88 - 96 and A = 110-116. Also, the 
SGII and Sk3 forces yield values of s comparable with 
each other that lie between the corresponding symmetry 
energy values when using SLy4 and LNS sets. The spe- 
cific nature of LNS force [48| (not being fitted to finite 
nuclei) leads to larger values of s (and to a larger size of 
the neutron-skin thickness, as it is seen from Fig. [2]) with 
respect to the results with other three forces. Although 
the values of s slightly vary within the Kr isotopic chain 
when using different Skyrme forces, the curves presented 
in Fig. [3] exhibit the same trend. 

The results shown in Fig. [3] are closely related 
to the evolution of the quadrupole parameter /3 = 
y^Q/(A(r 2 ) 1/2 ) (Q being the mass quadrupole mo- 
ment and (r 2 ) 1 / 2 the nucleus rms radius) as a function 
of the mass number A that is presented in Fig. |4] First, 
one can see from Fig. [H that the semi-magic A = 86 
and A = 118 Kr isotopes are spherical, while the open- 
shell Kr isotopes within this chain possess two equilib- 
rium shapes, oblate and prolate. In the case of open-shell 
isotopes, the oblate and prolate minima are very close in 
energy and the energy difference is always less than 1 
MeV. In this region of even-even Kr isotopes with very 
large N/Z ratio (> 1.7) the competition between the pro- 
late and oblate shapes has also been studied with HFB 
calculations and the Gogny force in Ref. [5lJ. Shape co- 
existence in lighter Kr isotopes has also been examined 
[52|,[53j]. Nevertheless, we specify in Fig. 2] which shape 
corresponds to the ground state of each isotope by en- 
circling them. Thus, the trend that the evolution of the 
symmetry energy shown in Fig. [3] follows can be clearly 
understood. The peaks of the symmetry energy corre- 
spond to the closed-shell nuclei that are spherical. Mid- 
shell nuclei (A = 96— 110) are well deformed and exhibit 
a stabilized behavior with small values of s. The tran- 
sitional regions from spherical to well deformed shapes 
correspond to transitions from the peaks to the valley in 
the symmetry energy. 

In Figs. [5] and [6] we give results for Sm isotopes 
(A = 140 — 156) as a well established example of de- 
formed nuclei. In the calculations, all Sm isotopes are 
found to have a prolate shape, except for the even-even 
144 Sm and 146 Sm nuclei that are spherical. Such an evo- 
lution of shape from the spherical to the axially deformed 
shapes in the same Sm isotopic region is in accordance 
with the results obtained from microscopic calculations 
in the RMF theory In Ref. [H| the ground state of 
the semi-magic 144 Sm (N — 82) is found to be spherical 
(having about 12 MeV stiff barrier against deformation) 
and the deformation in 146 Sm to be still small. With the 
increase of the neutron number, the ground state gradu- 
ally moves toward the deformed one till the well deformed 
i54-i58g m j^jj Also, the analysis of the potential energy 
curves [l5| within the same microscopic approach that 
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FIG. 1. (Color online) Proton (a) and neutron (b) rms radii differences relative to 86 Kr for oblate and prolate nuclear shapes, 
(c) Theoretical (with SLy4 Skyrme force) and experimental isotope shifts S < r\ > of Kr isotopes relative to 86 Kr. The results 
for oblate and prolate shape for A = 82, 84 isotopes are indistinguishable. 



we use in the present work, as well as within the HFB 
method with Gogny interaction (i~9j , confirmed the tran- 
sitional behavior between the spherical Sm and the 
well prolate deformed 154_158 Sm isotopes. 

The results for the correlation between the neutron 
skin thickness and the nuclear matter properties in fi- 
nite nuclei with SLy4 Skyrme force for a chain of Sm 
isotopes are shown in Fig. [51 while those with SGII, Sk3, 
and LNS forces are presented in Fig. [5] Similar to the 
case of Kr isotopes with transition at specific shell clo- 
sure, we observe a smooth growth of the symmetry en- 
ergy until the semi-magic nucleus 132 Sm [N — 82) and 
then an almost linear decrease of s while the neutron 
skin thickness of the isotopes increases. An approximate 
linear correlation between AR and po is also shown in 
Figs. [5jb) and[5Jb), while Fig. [5jc) exhibits a very ir- 
regular behavior of AR as a function of the asymmetric 
compressibility AK. Nevertheless, the values of AK de- 
duced from our calculations are in the interval between 
-295 and -315 MeV that compare fairly well with the 
neutron-asymmetry compressibility (K^ — —320 ± 180 
MeV) deduced from the data |55| on the breathing mode 



giant monopole resonances in the isotopic chains of Sm 
and Sn nuclei. 

The theoretical neutron skin thickness AR of Pb nu- 
clei (A = 202 — 214) against the parameters of interest, 
s, po, and AK , is illustrated in Fig. [7] In this work we 
consider an extended chain of Pb isotopes in compari- 
son to the one analyzed in Ref. [24| by adding two nuclei 
lighter than 206 Pb. Therefore, a more precise study of 
the corresponding correlations, especially in the transi- 
tion region at the double-magic 208 Pb nucleus, could be 
made. Here, we test three more parametrizatios (SGII, 
Sk3, and LNS) in addition to the SLy4 force applied in 
[24| . All predicted correlations manifest an almost lin- 
ear dependence and no pronounced kink at 208 Pb is ob- 
served. Similarly to Kr and Sm isotopes presented in this 
study (and isotopes from Ni and Sn chains described in 
[24j). the LNS force produces larger symmetry energies 
s than the other three forces also for Pb nuclei with val- 
ues exceeding 30 MeV. Another peculiarity of the results 
obtained with LNS is the almost constant AK observed 
in Fig. Et. 

Further attention deserves the value of the neutron skin 
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FIG. 2. (Color online) HF+BCS neutron skin thicknesses AR for Kr isotopes as a function of the symmetry energy s and the 
pressure po calculated with SLy4, SGII, Sk3, and LNS forces and for oblate and prolate shapes. The results for oblate and 
prolate shape for A = 82, 84 isotopes are indistinguishable. 



thickness in 208 Pb, whose determination has motivated 
recent experiments. The model-independent measure- 
ment of parity-violating asymmetry (which is sensitive to 
the neutron distribution) in the elastic scattering of po- 
larized electrons from 208 Pb at JLAB within the PREX 



Collaboration [27l . l28j has provided the first electroweak 
observation of Ai?=0.33t£;^ fm in 208 Pb. Obviously 
future precise measurements are needed to reduce the 
quoted uncertainties of AR. The distorted wave electron 
scattering calculations for 208 Pb [3l| extracted a result 
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FIG. 4. (Color online) The quadrupole parameter f3 as a 
function of the mass number A for the even-even Kr isotopes 
(A = 82- 120) in the case of SLy4 force. 



for the neutron skin thickness which agrees with that 
reported in the experimental paper (28{. As can be seen 
from Fig. Ufa), the value of AR for 208 Pb (0.1452 fm) de- 
duced from the present HF+BCS calculations with SLy4 
force agrees with the recent experimentally extracted skin 
thickness (0.1561^2? fm) using its correlation with the 
dipole polarizability (29J. However, this experimental 
value was derived by means of covariance analysis based 
on one Skyrme functional (SV-min). In this respect, a 
systematic study with a variety of EDFs as well as ex- 
perimental tests in other nuclei would be important be- 
cause the correlation between polarizability, neutron skin 
thickness, and symmetry energy is mo del- dependent (see, 



for example, Fig. 1 of Ref. [33[). In addition, our theo- 
retically obtained value of AR for 208 Pb agrees well with 
the value 0.18 ± 0.027 fm from Ref. Eel It is lower than 
the one obtained in Refs. [56[ and |57l| with the same 
Skyrme force, but is in agreement with the values cal- 
culated with self-consistent densities of several nuclear 
mean- field models (see Table I in Ref. [Hj]). The po and 
AK values for 208 Pb are in a good agreement with those 
from Ref. 

We would like to note that linear correlations AR vs. s 
and AR vs. po are found to exist in the Kr (A = 82 — 96), 
Sm (A = 140 - 156) and Pb (A = 202 - 214) isotopic 
chains. The correlation AR vs. AK presented on the 
examples of Sm and Pb isotopes (Figs. [5] and respec- 
tively) is less strong than the ones mentioned before and 
shows a slight irregular behavior. This concerns also the 
calculated results for AK in the case of Kr isotopes that 
are not shown in Fig. [2] Such observations are con- 
firmed by the results for these correlations obtained in 
Refs. [l^, In addition, it was demonstrated in [6lJ 
that among various observables that correlate with the 
neutron form factor (related to the neutron density and, 
thus, to the neutron skin thickness) the incompressibility 
is a poor indicator of isovector properties. 

Following our previous analysis within the CDFM ap- 
proach [24| . here we would like to give more detailed 
study of the weight function |/(a:)| 2 (that is related to 
the density and thus, to the structural peculiarities) to 
understand the kinks observed in the relationships be- 
tween AR and s, as well as AR and pq. The latter were 
shown to exist [24| in double-magic nuclei in the cases of 
Ni (at 78 Ni) and Sn (at 132 Sn) isotopic chains. As one 
can see in Figs. [2j [5j and|6]of the present work, they exist 
also in the considered cases of Kr (at 86 Kr) and Sm (at 
144 Sm) isotopes. In contrast, such a kink does not exist 
in the case of Pb isotopic chain (at 208 Pb, particularly). 

Here we analyze, as an example, the cases of Ni, Sn 
and Pb isotopic chains, trying to understand the origin 
of the kinks without additional complexities coming from 
deformation. For this purpose, let us introduce the quan- 
tity 



SA±2 - SA 
SA 



(30) 



which gives information on the relative deviation of the 
symmetry energy s of even-even isotopes with respect 
to the double-magic ones, namely with A = 78 for Ni, 
A = 132 for Sn, and A = 208 for Pb. Here we will con- 
sider the range of integration on x in Eq. (|10p in each of 
the cases of Ni, Sn, and Pb isotopes. Firstly, we will in- 
troduce the value of a; m j n at which the symmetry energy 
for nuclear matter s ANM {x) changes sign from negative 
(at x < x m i n ) to positive (at x > x m ; n ). In our work we 
use the symmetry energy s ANM (p(x)) [Eq. (fig)) ] from the 
Brueckner theory, that is defined by the second derivative 
of the energy per particle E(p 1 S). Considering in prin- 
ciple in the CDFM the range of x from zero to infinity 
we include in this way the region of densities po(x) from 



9 



0.16 
0.14 
0.12 
£ 0.10 
< 0.08 
0.06 
0.04 



b Sm 



SLy4 
(a) 



■ prolate 

■ spherical 



Sm 



u Sm 




L — prolate 
■— spherical 



bm 



27.6 28.0 28.4 28.! 

s [MeV] 



1.48 1.52 1.56 1.60 1.64 -330 -320 -310 -300 -290 -280 



p [MeV/fm 3 ] 



AK [MeV] 
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FIG. 6. (Color online) HF+BCS neutron skin thicknesses AR for Sm isotopes as a function of the symmetry energy s (a) and 
the pressure po (b) calculated with SGII, Sk3, and LNS forces. 



infinity to zero, respectively. When the values of x are 
small we consider in practice values of the density po(x) 
that are much larger than the density in the equilibrium 
state po- In this case unphysical (negative) values of the 
symmetry energy appear, and thus starting the integra- 
tion from x > x m i n we exclude these values. We have to 
note simultaneously that at x < x m i n the weight function 
|/(a:)| 2 is close to zero (it is its "left wing", see Fig. E}, so 
there is no contribution to s from this region (x < x m i n ). 
Secondly, we introduce the value of x max (in the "right 
wing" of |/(x)| 2 ) beyond which the contribution to s (i.e. 
the result of the integration in Eq. (|10p from x max to in- 
finity) is negligible. If we define by Ax = x max — £min, 



then we impose for x max the condition s — sax < 0.1 
MeV, where sax is obtained by Eq. (fl~0|) integrating over 
x from x m in to £ max , while s is the result of integration 
from x m - m to infinity. The points x m i n and £ max corre- 
sponding to the double-magic nuclei of the three isotopic 
chains considered are indicated in Fig. [5] where the in- 
trinsic weight functions \ f(x)\ 2 for 76 > 78 > 80 Ni [Fig.^a)], 

130,132,134 Sn [Fig.^b)], and 206,208,210p b [Fig. [^c)] , re- 

spectively, are separately presented. Then, if one uses 
Eq. (|10p that determines the symmetry energy s in finite 
nuclei within the CDFM, the contribution sax of s in the 
interval Ax is obtained. The values of sax are given in 
Table [J together with the values of cc m i n , x max , and s. We 
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can see from Table U that, as can be expected, the total 
symmetry energy is almost exhausted by its contribution 
sax, thus showing the important role which the nuclear 
surface plays for the obtained symmetry energy values. 
For a better illustration of this fact, we plot together in 
Fig.Hthe weight function \f(x)\ 2 and the HF+BCS total 
density p (R) of 208 Pb in the case of SLy4 force. As can 
be seen, the nuclear surface part of the density which de- 
termines the weight function |/(x)| 2 containing the peak 
around the maximum is responsible to a large extent for 
the main contribution to the total symmetry energy. 

The analysis of the quantity As± [Eq. (|30|) ] might be 
instructive because this quantity is a direct measure of 
the relative deviation of the symmetry energy with re- 
spect to the double-magic nuclei taking them as reference 
nuclei in each of the chains, where the kinks are expected. 
The values of As + and As_ are listed in Table [TTT where 
the two numbers for each isotopic chain correspond to 
the range of integration Ax. One can see firstly from 
this Table that the absolute values of As+ and As_ for 
Pb isotopes are comparable with each other, which is not 
the case for the two other isotopic chains. Secondly, and 
very important is that the As + value turns out to be 
negative and As_ value to be positive for Pb isotopes 
at the range of integration Ax, and this is the main dif- 
ference regarding to the corresponding values (both are 
negative) in the Ni and Sn chains. 

These differences can be attributed to the profiles of 
the density distributions, particularly in the surface re- 
gion. They are given in Fig. 1 1 01 where curves for five Ni, 
Sn, Pb isotopes around double-magic 78 Ni, 132 Sn, and 
208 Pb nuclei are presented in panels (a), (b), and (c), re- 
spectively. One can see from Fig. [TU] the same trend in 
the tails of the three isotopic chains, which are ordered 
according to the mass number A, being higher for heav- 
ier isotopes to produce larger radii. On the other hand 
the behavior in the top part of the surface region, shown 



in the inset of the panels, is different. In the case of Ni 
and Sn isotopes in panels (a) and (b), one observes that 
the double-magic nuclei have the largest density with all 
the neighboring isotopes lying below. In the case of Pb 
isotopes in panel (c), the density increases from heavier 
to lighter isotopes with the double-magic nucleus in be- 
tween. In Pb isotopes, this ordering is opposite in the 
tail. As a result of this, the slope of the density in Pb 
isotopes, and therefore |/(x)| 2 , decreases with the num- 
ber of neutrons continuously and no kink is present in the 
symmetry energy. On the other hand, in the case of Ni 
and Sn isotopes, the slope of the density is larger for the 
double-magic isotopes generating a kink in the symmetry 
energy. 



In this way the kinks displayed in our previous study 
[24j by the Ni and Sn isotopes and the fact that no kinks 
appear in the Pb chain considered can be understood. 
It concerns not only the symmetry energy evolution, but 
also its relationship with the neutron skin thickness AR 
and the pressure pq. We would like to summarize the 
discussion about the kinks, saying that within the CDFM 
we are able to use only the densities resulting from the 
specific shells that are occupied. The presented Figures 

and ITD1 give information about the differences between 
densities of Ni, Sn, and Pb isotopes that are reflected 
in the corresponding weight functions |/(x)| 2 for these 
isotopes (see Fig. [5]). Of course, these differences are due 
to the different occupancies of the shells in these nuclei. 
Thus, we see that in general kinks are produced at shell 
closures, but the analyses of the precise dependence of 
various kinks on the amount of occupation of specific 
shells will require further work. 
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130, 132, 134) isotopes (b), and three Pb (A = 206, 208, 210) isotopes (c) by using the HF+BCS total densities for these nuclei 
and with SLy4 force. 



TABLE I. Integration limits values s m i n and i max (in fm), values of the contribution sax to the total symmetry energy s (in 
MeV) for three Ni (A = 76, 78, 80), Sn (A = 130, 132, 134), and Pb (A = 206, 208, 210) isotopes. 





76 Ni 


78 Ni 


80 Ni 


130 Sn 


132 Sn 


134 Sn 


206p b 


208 pb 


210p b 


^min 


3.26 


3.30 
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X max 
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FIG. 9. The weight function \f(x)\ 2 and the HF+BCS total 
density po(R) of 208 Pb calculated with SLy4 force. 



TABLE II. Relative deviation values of the symmetry energy 
As+ and As_ [Eq. P0[) ] for the range of integration Ax in 
Eq. (|10[) and for Ni, Sn, and Pb isotopes. 





Ni 


Sn 


Pb 




-0.0137 


-0.0070 


-0.0035 


As_ 


-0.0072 


-0.0049 


0.0038 



IV. CONCLUSIONS 

In this paper, we have investigated possible relation- 
ships between the neutron skin thickness of deformed 
neutron-rich nuclei and the symmetry energy character- 
istics of nuclear matter for these nuclei. A microscopic 
approach based on deformed HF+BCS calculations with 
Skyrme forces has been used. Four Skyrme parametriza- 
tions were involved in the calculations: SGII, Sk3, SLy4, 
and LNS. Nuclear matter properties of nuclei from Kr 
and Sm isotopic chains have been studied by applying 
the CDFM that provides a transparent and analytic way 
to calculate the intrinsic EOS quantities by means of a 
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FIG. 10. (Color online) HF+BCS total densities in the surface region for five Ni (A=74-82) isotopes (a), five Sn (A=128-136) 
isotopes (b), and five Pb (A=204-212) isotopes (c) around double-magic 78 Ni, 132 Sn, and 20S Pb nuclei, respectively, calculated 
with SLy4 force. 



convenient approach to the weight function. As a first 
step we study the variation with neutron number of pro- 
ton and neutron radii predicted by the self-consistent mi- 
croscopic (HF+BCS) calculations. We find that charge 
and neutron radii increase similarly with neutron num- 
ber following the increase of neutron radii, as it should be 
expected from the general properties of the nuclear force. 
The same microscopic method is then applied to calculate 
the skin thickness and the weight function for each iso- 
tope. The analysis of the nuclear symmetry energy s, the 
neutron pressure po, and the asymmetric compressibility 
AK has been carried out on the basis of the Brueckner 
EDF for infinite nuclear matter. 

For both Kr (A = 82 - 96) and Sm (A = 140 - 156) 
isotopic chains we have found that there exists an approx- 
imate linear correlation between the neutron skin thick- 
ness of these nuclei and their nuclear symmetry energies. 
Comparing with the spherical case of Ni, Sn, and Pb 
nuclei described in our previous study 241 ]. we note that 
the linear correlation observed in the Kr and Sm isotopes 
is not smooth enough due to their different equilibrium 
shapes, as well as to the transition regions between them. 
As known, the latter are difficult to be interpreted as they 
exhibit a complicated interplay of competing degrees of 
freedom. Nevertheless, a smoother behavior is observed 
in Kr isotopes that is a consequence of the stabilization 



of the oblate shapes along the isotopic chain. As far 
as Sm isotopes are concerned, the shape evolution from 
the spherical to the axially deformed configurations in 
the Sm isotopes causes a less pronounced linearity of the 
observed correlation between AR and s. A similar cor- 
relation between AR and po is also found to exist, while 
the relation between AR and AK exhibits an irregular 
behavior. However, for both classes of deformed nuclei 
an inflection point transition at specific shell closure, in 
particular at semi-magic 86 Kr and 144 Sm nuclei, appears 
for these correlations of the neutron skins with s and po- 
In addition, the role of the relative neutron-proton asym- 
metry on the evolution of the symmetry energy has been 
pointed out on the example of Kr isotopes with N > 60. 

We have analyzed in detail the existence of kinks on 
the example of the Ni and Sn isotopic chains and the 
lack of such kink for the Pb isotopic chain. For this pur- 
pose, we have estimated the relative deviation of the sym- 
metry energy of even-even isotopes with respect to the 
double-magic 78 Ni, 132 Sn, and 208 Pb nuclei. An anal- 
ysis of the weight function |/(a;)| 2 (determined by the 
HF+BCS density distributions) as the key ingredient of 
the CDFM is performed. It is shown that for the Pb 
isotopes the different signs of the relative deviations cor- 
responding to the range of integration on x in Eq. (|10|) 
that contains the peak of |/(a;)| 2 is in favor of the absence 
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of kink in the Pb chain. Thus, from the previous study 
in Ref. [24j and the present analysis the kinks displayed 
by the Ni and Sn can be understood as consequences of 
particular differences in the structure of these nuclei and 
the resulting densities and weight functions. 

It has to be mentioned that the used microscopic theo- 
retical approach is capable also to predict important nu- 
clear matter quantities in deformed neutron-rich exotic 
nuclei and their relation to surface properties of these nu- 
clei. This is confirmed by the good agreement achieved 
with other theoretical predictions and some experimen- 
tally extracted ground-state properties. New experimen- 
tal results on giant resonances, neutron skin in heavy 



nuclei, and heavy-ion collisions could further lead to new 
constraints, e.g. on the nuclear symmetry energy, allow- 
ing successful interpretation of data and observations on 
physical quantities of nuclear systems. 
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